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ABSTRACT 

Tidal dissipation may be important for the internal evolution as well as the 
orbits of short-period massive planets — hot Jupiters. We revisit a mechanism 
proposed by Ogilvie and Lin for tidal forcing of inertial waves, which are short- 
wavelength, low-frequency disturbances restored primarily by Coriolis rather than 
buoyancy forces. This mechanism is of particular interest for hot Jupiters be- 
cause it relies upon a rocky core, and because these bodies are otherwise largely 
convective. Compared to waves excited at the base of the stratified, externally 
heated atmosphere, waves excited at the core are more likely to deposit heat in 
the convective region and thereby affect the planetary radius. However, Ogilvie 
and Lin's results were numerical, and the manner of the wave excitation was 
not clear. Using WKB methods, we demonstrate the production of short waves 
by scattering of the equilibrium tide off the core at critical latitudes. The tidal 
dissipation rate associated with these waves scales as the fifth power of the core 
radius, and the implied tidal Q is of order ten million for nominal values of the 
planet's mass, radius, orbital period, and core size. We comment upon an alter- 
native proposal by Wu for exciting inertial waves in an unstratified fluid body 
by means of compressibility rather than a core. We also find that even a core of 
rock is unlikely to be rigid. But Ogilvie and Lin's mechanism should still operate 
if the core is substantially denser than its immediate surroundings. 

Subject headings: hydrodynamics — waves — binaries: close — planetary systems — 
stars: oscillations, rotation 



1. Introduction 



The discovery of extrasolar planets has sharpened the need for a predictive theory 
of tidal circularization and synchronization. Some 2% of nearb y single FGK stars ha rbor 
roughly Jupiter-mass planets with orbital periods below ten days (jCumming et al.ll2008l . and 
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references therein). The orbital eccentricities of the shorter-period planets, especially those 
below four days, are mark edly less than thos e of the longer period systems, presumably as a 
result of tidal dissipation (IRasio et al.lll996l ). Synchronization is likely to occur more easily 
and at longer periods than circularization because of the sma ller moment of inertia associated 
with the spin as compared to the orbit (ILubow et al.lll997l ). For close stellar binaries, the 
lack of an adequate tidal theory causes little uncertainty concerning the internal evolution 
of the components since tidal heating contributes negligibly to the luminosity, and since it 
is assumed that tides are always sufficient to synchronize the spins of stars that come into 
contact with their Roche lobes. The intrinsic luminosities of jovian planets due to their 
gradual contraction and loss of primordial heat are so small, however, that tidal heating 
may be competitive. Indeed tides have been invoked to account for the anomalously large 
radii and low densities, compared to baseline mo dels, of some planets that are observed to 
transit their host stars (IBodenheimer et al.ll200ll ). 



Most studies of orbital evolution reduce the tidal uncertainties to a single parameter, 
the tidal quality factor Q, which is inversely proportional to the tidal dissipation rate, and 
which one hopes to calibrate by reference to stellar binaries of known age, such as those in 
star clusters, or to the inferred tidal in teractions between Jupiter and its Galilean satellites 
( IGoldreich &: Soterl Il966l ; iMazehl 120081 . and references therein) . A difficulty with such pure 
empiricism is that Q may depend upon structural details such as composition, equation of 
state, rate of rotation relative to the body's density or to the tidal period, and so on, which 
differ between the object of interest and the rather limited set of calibra tors. Indeed there 
is some indirect evidence that Q does vary. As has been pointed out (IRasio et al.l Il996 



Sasselovl 120031 ; IJackson et al.l 120081 ) , differences in Q between planets and stars, as well a s 
possible dependence of Q on the ratio of the orbital to spin frequency (jOgilvie fc Linll2007l ). 
may be crucial to the s urvival of hot Jupiters: since the host stars rotate subsynchronously 
( iFabrycky et al.ll2007al ) and the mass ratio is large, tidal dissipation within the star would 
tend to drag the planet inward. 

In order to assess the importance of tidal inputs for hot Jupiters, furthermore, it is 
necessary to predict not only the overall rate of tidal dissipation, but also where within 
the planet mechanical energy is converted to heat. This requires additional assumptions 
or knowledge beyond an empirically calibrated Q alone. For example, perhaps the best 
understood and most predictive tidal mechanism, at least for nonrotating bodies, is the 
excitation of shortwavelength g-modes at an interface between convective (isentropic) and 
radiative (stratifie d) regions, originally proposed by Zah n for applicat i on to early- type stars 
(Zahn 197ol . 1 19751 ) and first applied to hot Jupiters by Lubow et al. ( 1997 ). The g-modes 
would dissipate within the radiative region since they do not propagate in isentropic regions. 
For irradiated planets, as for early-type stars, this means that dissipation would occur in the 
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outer parts. On the other hand, if turbulent viscosity due to convection were effective (IZahn 
19661 ). then the heat would be deposited locally at depth. Because the thermal timescales 
and density scale heights differ greatly between the convective and radiative regions, the g- 
mode and turbulent mechanisms would have different consequences for the planetary radius 
even if they produced the same Q. 



As argued by lGoldreich fc Nicholson! (119771 ). the effective viscosity of convection is prob- 
ably very strongly suppre ssed when the turnover time of the larger convective eddies exceeds 
the tidal period (see also lZahnlll989l ; iGoodman fc Ohlll997l ; iPenev et al.ll2007l ). Since jovian 
planets are very deeply into this regime (r conv ~ 10 3 yr), it is unlikely that convection can be 
responsible for Q values as low as are inferred for Jupiter, Q p « 10 5 — 10 6 (IGoldreich fc Soter 
19661 ; iPeale fc Greenberdll980l ). While bearing in mind that phase transitions or other non- 



trivial microphysics may contribute to tidal dissipation (e.g. iStevensonl ll98Cj ) . we take the 
traditional view that the likely alternative is a dynamical tide: that is, resonant tidal forcing 
of a low-frequency, short-wavelength mode or traveling wave at special locations within the 
planet, such as the convective-radiative interface already discussed. 

Short waves have two basic physical advantages as candidates for a tidal theory. First, 
they are more easily d amped tha n the large-scale ellipsoidal dis tortion ( "equilibrium tide" ) , 



by radiative diffusion (jZahnlll975l ) or by nonlinear breaking (e.g. IGoodman fc Dicksonlll998l ). 
Second, the number of short-wavelength modes is potentially large, scaling as (R p /X) 2 for 
modes with a fixed azimuthal dependence exp(im<f)). For direct forcing, the frequency of 
the mode as well as the azimuthal order m must match those of the tidal potential. Viewed 
in the corotating frame of the planet (we assume a uniformly rotating background state), 
this frequency is normally much lower than the fundamental dynamical freqency Wdyn = 
(3GM p / Rp) 1 / 2 , so the modes of interest must be approximately noncompressive and restored 
by buoyancy and/or rotation rather than pressure. The nomenclature for such low-frequency 
motions is rich and bewildering: internal waves, g modes, r modes, toroidal modes, hybrid 
modes, Hough modes. These names distinguish the relative importance of buoyancy versus 
rotation and other special properties. In this paper, we concentrate on modes restored by 
rotation rather than buoyancy, which we refer to collectively as "inertial" modes or waves. 

All of these nearly incompressible motions have frequencies in limited ranges controlled 
by the Coriolis and Brunt- Vaisala parameters, regardless of wavelength. Local dispersion re- 
lations for the wave frequency depend upon the ratios of components of the three-dimensional 
wavevector, but hardly depend upon the wavelength itself, except via viscous or radiative 
damping terms. This behavior is completely different from that of modes depending upon 
compressibility (acoustic or p modes) or selfgravity. In the ideal-fluid limit, the spectrum 
of modes that can be resonant with a tide in the allowed frequency range is therefore dense 
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( IPapaloizou fc Pringldll98ll ). But, the tidal potential must have a finite projection onto a 
resonant mode in order to excite it. So, the main challenge for a theory of dynamical tides 
is to estimate these projections, also called overlap integrals. 

Two recent and independent studies have proposed that inertial modes may be tidally 
forced even in the absence of a stably stratified surface layer, and at rates that approach 
what is needed to explain the observationally inferred Q p , provided that the tidal frequency 
in the corotating frame is less than twice the rotation frequency; this is the relevant regime 
for circulariza t ion o f a planet that is already synchronized. By direct numerical methods, 
Ogilvie &: Linl (120041 . hereafter OL04) calculated the linear excitation of inertial waves/modes 
in a compressible spherical annulus surrounding a solid core. The current belief is that 
Jupiter itself contains 15 — 30M e of heavy elements (atomic number Z > 2), or 5 — 10% of 
the planet's total mass, and also that many or most of the hot J up iters may be ev en more 



enriched; whether these "metals" are concen trated in a distinct core is uncertain (jGuillot 



20051 ; iGuillot et al.ll2006l ; iBurrows et al.ll2007l ). OL04 calculate the tidal dissipation in steady 
state by balancing the excitation against an artificial viscous term, but they conclude that 
the tidal dissipation due to the excitation of inertial modes, though varying erratically with 
frequency in the inertial range (i.e., |w t ide I < |2fi|), tends towards finite values corresponding 
to Q ~ 10 6 in the inviscid limit. This they explain by the presence in their low- viscosity 
models of very short-wavelength disturbances concentrated on a "web of rays" in the poloidal 
plane. They find that the intensity of this short-wavelength response correlates with the size 
of the core, and they suggest that it has something to do with closed ray paths ("wave 
attractors" ) reflected alternately by the inner and outer boundaries of the spherical annulus. 
OL04's i nterpretation of their own num erical r esults seems to have be en informed by previous 
work by lRieutord fc Valdettarol (119971 ) and by lRieutord et al.l (120011 ). who studied the linear 
modes of an incompressible, slightly viscous liquid in a spherical shell, by a combination 
of numerical and analytic methods. The latter authors demonstrated the existence of very 
fine features in the velocity, and they explored the relationship of these features to critical 
latitudes and wave attractors, but they did not consider tidal excitation. Subsequently, 
Ogilvid (120051 ) used a simplified analytical model to argue that wave attractors can absorb 
energy at nonzero rates that vary continuously with forcing frequency in the inviscid limit. 

The models of OL04 include a stably stratified atmosphere, with an apparently inde- 
pendent set of "Hough modes" excited at the interface between this zone and the convective 
region, as previously demonstrated for rotating massive s tars by Savonije and Papaloizou 
(ISavonije &: Papaloizoulll9971 ; iPapaloizou fc Savonijdll997l ). The Hough modes appear to be 
much less dependent on the core, and the dissipation associated with them varies much more 
smoothly with tidal frequency; also, these modes extend outside the inertial range, i.e. to 

Ntide > |2f2|. 
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Shortly after the work of OL04, Y. Wu claimed to de monstrate th e tidal excitation of 
inertial mode s in enti r ely un stratified and coreless bodies. In lWul ( I2005al ). exploiting methods 
developed by iBryanl (118891 ) for incompressible rotating bodies, she analyzes the properties 
of free (unf orced) mode s of oscillation in compressible models with special radial density 
profiles. In IWul (j2005bl ). she calculates spatial overlap integrals between these modes and 
a quadrupolar perturbing potential to find Q ~ 10 9 ; she argues that Q may fall to < 10 7 
in more realistic models with a radial density jump due, for example, to a first-order phase 
transition. Usi ng a different mathematical formalism, but similar underlying low-frequency 
approximation, llvanov &: Papaloizoul (120071 ) have calculated tidal excitation of inertial modes 
in planets on highly eccentric orbits, with realistic but isentropic equa tions of state. This is 
an ext ension of earlier calculations for rotating n = 3/2 polytropes by iPapaloizou fc Ivanov 
( 120051 ). Their models are coreless, like those of Wu, but unlike hers, their tidal excitation 
is dominated by two large-scale inertial modes, perhaps because they treat the perturbing 
potential as uncorrelated from one periastron encounter to the next, so that the excitation 
is non-resonant. 

It is clear from the above that considerable progress has been made in recent years 
on the dynamical tides of rotating planets, but the importance of short-wavelength inertial 
waves remains obscure: in particular, the circumstances under which they can be tidally 
excited. OL04's results indicate that a solid core is somehow important, while Wu's work 
suggests that it may not be. It is desirable to achieve a semianalytic understanding of the 
excitation before making detailed calculations for realistic planetary interiors. For one thing, 
since the structure of hot Jupiters is still much less well understood than that of stars, it will 
be useful to know what features of the structure are most important for the tidal problem. 
For another, the apparently chaotic variation of tidal torqu e with tidal frequency seen i n 
the results of OL04, and to some extent in the earlier ones of ISavonije fc Papaloizoul (119971 ). 
raise questions about what is required of a numerical calculation to obtain convergence in 
the inviscid limit, or indeed what it means to converge. 



2. Basic Equations 

We assume an isentropic body in uniform rotation at angular velocity ft. The linearized 
equations of motion in the corotating frame are 

-iuv + 2SlXv = -V (c s 2 — + $! ] = -VV>, (T 



P 

-%ujp x + V • (pv) = 0. (2) 
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Where necessary, first-order eulerian perturbations have been distinguished by a subscript 
from the corresponding quantities in the background state: for example, first-order mass 
density p\ and pressure p\. All such perturbations have the time dependence exp(— iut), 
uj being the angular frequency of the tide viewed from the corotating frame. Since the 
unperturbed fluid velocity vanishes in this frame, v is understood to be of first order even 
though it is not explicitly marked as such. We have introduced %p as the sum of the first-order 
enthalpy Pi/p = c^pi/p and gravitational potential $1, where c 2 = (dp/dp) 5 is square of the 
sound speed. 

The gravitational potential can be further subdivided as $ = $i, ex t + $i,seif, the first 
term representing the tidal potential exerted by the companion, and the second representing 
the perturbation to the self-gravity of the fluid body. For a fully self-consistent treatment, 
the equations above should be supplemented by Poisson's equation V 2 $i, S eif = AnGpi, but it 
will not be necessary to deal with $i jSe if explicitly in the present paper, since we are mainly 
concerned with short- wavelength inertial modes for which self-gravity is unimportant. Even 
when we focus on the large-scale ellipsoidal distortion of the body, for which the self-potential 
is important, we pretend that it has already been computed and included in $1. 

Equation ([1]) can be solved algebraically for the velocity: 

v = — — — - + 20 x + (iu)-Hnn • wl = m w • V^. (3) 

A\V — cu z 

This defines M w as 3 x 3 matrix or tensor that is spatially constant in cartesian coordinates. 
Using ([3]) to eliminate v from the linearized continuity equation (|2J), and remembering that 
tp = $1 + CgPx/p leads to a wave equation for ip: 



4fi 2 -^ 2 _ 1 
(ip-$i) = -V 

C s P 



2 40fi 
p ( + — O x Vip — • 



(4) 



For application to hot Jupiters and binary stars, the tidal frequency uj is typically comparable 
to the rotation frequency Q, and both of these are small compared to the dynamical frequency 
co>d y n = (3GM p /_Rp) 1 / 2 . Throughout most of the planet, c s ~ uJd yn R p . It follows that the 
lefthand side of eq. (j3J) is negligible throughout most of the interior for disturbances whose 
wavelength A is small compared to the planetary radius R p . However, the lefthand side is 
important near the surface of the planet; in fact it is singular there, because c 2 — > at the 
surface, at least with idealized zero-temperature, zero-pressure boundary conditions. The 
righthand side of eq. (jl]) contains a term oc v • V In p that is also singular at the surface. 
In order that there be a well-behaved solution for ip and v, it is necessary that these two 
singular terms should balance: v • V lnp — > iuj(ip — $i)/c 2 = iujp\j p. Using the hydrostatic 
equilibrium of the unperturbed state, this can be recast as 

iuj 

v • n = — ($1 — ip) at the surface, R = R p , (5) 
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where n is the outward-pointing normal to the unperturbed boundary and g > is the 
effective gravity. This is the usual free boundary condition: it says that the lagrangian, 
not eulerian, enthalpy perturbation vanishes at the surface. The normal components of 
the velocity and of the displacement £ = v/(—iuj) do not vanish at the boundary. It is 
true that it is possible to divide the tidal response into long-wave and short-wave parts, 
v = v\ ong + fghort, in such a way that it would be an excellent approximation to neglect 
n * v short at the boundary. But then the analog of equation for -^short would contain 
inhomogenous terms involving ^long even if the terms oc c~ 2 were neglected. While the 
principle of the free boundary condition is familiar, w e have emp hasized the point because 



it will be important to our discussion of the results of IWul (l2005bl ) in §4.11 

For the time being, we represent the rocky core by a rigid sphere of radius R c < i? p : 

v ■ h = at R = i? c , (6) 

where n is now the normal to the core. We take n to point away from the fluid, that is, 
downward at the core and upward at the surface. In Q4.'2\ we show that that the core should 
deform with the equilibrium tide, and that it is more likely to be liquid rather than solid. 
To accomodate this deformation, the tidal Q derived for a completely rigid core [eq. (138]) ] 
requires an overall correction factor of order unity that depends upon the density contrast 
between the core and the convective region. 



2.1. Incompressible limit 

Many of the calculations of this paper will be carried out in the limit c 2 — > oo. For a 
consistent hydrostatic equilibrium in the unperturbed state, the background density must 
be constant. In this limit, eq. (j3j) simplifies to 

V 2 ^-(2fi- V) 2 ^ = 0, 

d 2 ip d 2 ip 1 d 2 tp 
dx 2 dy 2 c 2 dz 2 

in cartesian coordinates with the z axis parallel to CI. We have introduced the abbreviation 

- U (8) 
yjAVL 2 - U 2 

for the coefficient of the vertical derivatives. It is dimensionless, unlike the sound speed 
c s . When u 2 < 4f2 2 , equation jTJ is hyperbolic, z plays the role of a timelike coordinate, 
and c clearly plays the role of wave "speed." The hyperbolic nature of the more general 
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equation (0J in the low-frequency regime has been noted by ISavonije fc Papaloizoul (119971 ) 
and emphasized by OL04. 



Equation (j7j) has been studied extensively for rotating incompressible fluids (e.g. [Greenspan 



19691 ). In the incompressible limit, the potential perturbation does not enter the wave equa- 
tion at all. But it does enter the free boundary condition (151) . which is still applicable. 

A constant-density body is not a realistic model for a planet or star. However, in 
addition to simplifying calculations, this model exhibits particularly clearly the division of 
the tidal response between short and long-wavelength parts, which is the main object of this 
paper. 



2.2. WKB dispersion relation and group velocity 

For waves sufficiently short that they may be described locally by plane waves propor- 
tional to exp(ik ■ x — iuit), and for frequencies uj 2 <C c 2 k 2 , equations PJ or ([7]) lead to the 
dispersion relation 

where k = \k\. For free oscillations, presuming that the components of the wavevector k are 
at least approximately real (so that the envelope of the wave varies slowly with position), 
the dispersion relation requires that — 2Q < uj < 2Q. It is important that uj is independent 
of wavelength A = k/2ir; it depends only on the direction k = k/k, which must lie on a 
double-napped cone whose axis is vertical, i.e. parallel to fi, and whose half angle is 



(10) 



Hereafter, to fix important signs, we adopt the conventions uj > and Q > 0. Tidal 
components that are retrograde with respect to the planetary spin will be represented by 
negative values of the azimuthal quantum number m rather than negative values of uj: that 
is, nonaxisymmetric tides have the ((f), t) dependence exp(zm0 — iujt) with and their 

azimuthal pattern speed is uj/m. With this convention, uj = \2fl • k\/k, and the group 
velocity becomes 

duj . k x (2Q x fc) 

gv= dk = Slgn( z) k • ( } 

So the waves and the energies they carry move crabwise, at right angles to their wavevectors. 
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The time-averaged energy density and energy flux carried by inertial waves are 

£ - i^- = 2(ifei*i' (12a) 

T = £V gT) (12b) 
if the physical velocity perturbation is the real part of v. 



2.3. Tidal response of a coreless incompressible planet 

Inasmuch as the planet is small compared to its orbital semimajor axis (a) — typically 
R p /a ~ 10 -2 for hot Jupiters — the dominant component of the external tidal potential is 
quadrupolar, that is, proportional to R 2 Y 2)m (9, <f>) exp(-iut) in spherical polar coordinates 
R9(f> centered on the planet. Expressed in cartesians, $i )ex t is then a linear superposition of 
the terms 

m = ±2: A± 2 (x±iy) 2 (13a) 
m = ±1 : A±iz(x±iy) (13b) 
m = 0: A (2z 2 - x 2 -y 2 ), (13c) 

times exp(— iut), where the coefficients A± 2 , A±\ and Aq are complex constants. On the 
other hand, it is clear that the wave equation (j3J) can be satisfied by polynomials in the 
cartesian coordinates. If the unperturbed body is axisymmetric — as it should be when 
Q 2 <C uj 2 yn ) — then after transients have died away, the forced responce of if) must have 
the same azimuthal symmetry as that of $i )ex t- If this response is also of second degree 
in the coordinates (an assumption that turns out to yield an acceptable solution), then 
its components at m = ±1 and at m = ±2 must have the same form as the first two of 
equations ({TBI ; it is easily seen that these satisfy the wave equation (J7J). But the axisymmetric 
component must have the form 

m = : ip oc 2c z + x + y + constant, (14) 

which differs from that of the perturbing potential, eq. (113c| . 

Following eq. ([31 , the displacements associated with these quadratic forms of if) are linear 
functions of the coordinates. From this and the constancy of the density in the interior, it is 
easily shown that $i )Se if oc $i, ex t) so that the total potential perturbation $i = $i )Cxt + $i, S eif 
is also a superposition of the terms (fTBT) . with an appropriate rescaling of the coefficients. 

It remains to relate the amplitudes of if) and $i using the surface boundary condition 
(JSJ). If we neglect the influence of centrifugal force on the unperturbed state on the grounds 
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that Q 2 <C u djn = 4:7iGp, the unperturbed boundary is spherical, x 2 + y 2 + z 2 = R 2 , so that, 
with use of eq. (j3J), 



v • n = v r 



2iQc 2 

zd z - c 2 (xd x + yd y ) H (yd x - xd y ) 



ip- (15) 



iuR 

The solutions for the azimuthal harmonics of ip then turn out to be 

* - s^S*- «»^ <-> 

* - g^ff^-lo,) * ifme±1 ' (16b) 

^ = ^^ + ^ + 2cV 1 ifm = 0. 16c 

These have been calculated to lowest order in u 2 /nGp and Q 2 /7iGp, which are much less 
than unity. At this level of approximation, ip is negligible on the righthand side of the surface 
boundary condition ([S]), which therefore reduces to eq. (fT7|) . exactly as for the traditional 
equilibrium tide when rotation is neglected. Also to this order in Q 2 /7iGp, the boundary is 
spherical. A constant term in ip is not required even when m — 0, since the radial velocity and 
radial displacement derived by substituting eq. <HM into eq. (TTBT) turn out to be oc P 2 (cos9); 
this is no miracle, but rather a consequence of the fact that 

^«4x [l + 0(Q 2 /nGp)}, (17) 
9 

where g is the angular average of the surface gravity (= 4irGpRp/3 for a const ant- density 
planet). 

Inasmuch as u and therefore c are real, the responses (FIB]) are perfectly in phase with 
the tidal forcing, so that there is no secular input of tidal energy to this (highly idealized) 
planet. An exception might occur if the denominator in eq. ffl6l) were to vanish, when the 
tide would be in resonance with an (£,m) = (2, 1) free-precession mode. But this happens 
only at the single discrete frequency uj = Q [in inertial space, u = 20] , and then only if the 
spin and orbit are misaligned. So it cannot be a general explanation for tidal dissipation. 
Also, the responses ffTBl) are entirely smooth and long-wavelength — or rather, nonwavelike, 
since the individual components of ip in (fl6|) have no radial nodes at R > 0, as is the case 
for $i itself. Short-wavelength inertial oscillations are not excited. Note that the inner 
boundary condition (jSj) is satisfied because we are assuming that R c = 0, and because all 
components of the velocities, being linear functions of the coordinates, vanish at the origin. 

These results are hardly new. Solutions f or the circular but not necessarily synchronized 



case are known as Roche- Riemann ellipsoids (I Chandr asekharl 1 1 98 71 ) . Viewed in the corotat- 



ing frame of the orbit rather than that of the body's spin, they appear as ellipsoids with 
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stationary axes, though the flow velocity is generally nonzero in this frame. The axis ratios 
of these classical solutions are not limited to values near unity, as ours are; they are e xact 
nonl inear solutions. Though ex act, they are not always stable ( iChandrasekharl 119871 ). In 
fact iLebovitz fc Lifschita (119961 ) have shown that except for the synchronous solutions, the 
Riemann and Roche-Riemann ellipsoids are generically vulnerable to small-scale parametric 
instabilities, with growth rates that increase with increasing departure from stationary and 
circular streamlines. However, the growth of small-scale inertial waves by parametric insta- 
bility is entirely distinct from direct tidal forcing, and the secular energy dissipation rate that 
results is much more difficult to estimate because it depends upon the amplitu des at which 



these instabilities saturate, which necessarily involves nonlinear considerations (lArras et al. 
20031 . and references therein). 



We conjecture that the absence of small-scale waves from the tidal response is not a 
peculiarity of the incompressible, constant-density limit, but that it holds more generally for 
isentropic, coreless, compressible bodies whose unperturbed enthalpy profiles are sufficiently 
smooth, provided that t he outer bo undary is free. This conjecture, if true, would seem to 
contradict the results of IWul (]2005bl ). More will be said on this subject in §4.11 



3. Production of short waves by scattering from the core 

When the radius of the rigid core is nonzero, it does not seem to be possible to con- 
struct a purely nonwavelike tidal response, confirming what OL04 concluded from numerical 
computations with a compressible model. As far as we can tell, this remains true even in our 
incompressible, const ant- density model, at least in the limit that the tide is strictly periodic. 
It seems that it ought to be possible to prove (or disprove) this rigorously We have no such 
proof, but we explain briefly why we think it is true in §3 .It the material of that subsection 
is not used in the rest of the paper. Then in §3.2[ we go on to the main goal of this section, 
which is to estimate the production of short waves at the core on the assumption that the 
response is partly wavelike. 



3.1. Homogeneous solutions of the wave equation 

The difficulty in constructing a purely nonwavelike response appears to lie in the lack of 
a sufficiently rich set of homogenous functions that solve equation (j4j) and are nonsingular 
on the surface of a sphere, or indeed any closed surface. A function f(x,y,z) is said to be 
homogenous of degree D if / (ax, cry, a z) = cr D f(x,y, z) for any constant a ^ 0; the degree 
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may be negative. A homogeneous function cannot have isolated radial nodes, since if it 
vanishes at any point, then it must vanish at all points along the same radial ray. But of 
course a linear combination of such functions may have radial nodes. 

If c 2 < 0, as when u 2 > AQ 2 [cf. eq. ([8], then equation (J4]) is elliptic rather than hy- 
perbolic, and after the rescaling z = z/y/—c 2 , reduces to Laplace's equation. Solid spherical 
harmonics r e Ye tm (9, 4>), when expressed in cartesians, are homogeneous solutions — in fact 
polynomials — of degree D = i > 0. Each of these functions has a companion r^^ 1 !^^, 0) 
with negative degree, D — — (£+ 1), that also solves Laplace's equation. Given smooth val- 
ues of ip or dip/dn on concentric spheres or (less conveniently) on spheroids, it is possible to 
construct a smooth solution of Laplace's equation in the space between them as linear com- 
bination of these functions; in general, unless the boundary conditions are very restricted, 
both signs of D are required. 

As is well known, the spherical harmonics of a given degree I constitute an irreducible 
representation of the rotation group SO(3), to which the laplacian is invariant. This suggests 
that when uo 2 < 4f2 2 , we should be concerned with homogeneous functions that belong to 
irreducible representations of SO(2, 1), the 2+1- dimensional version of the Lorentz group, 
since that is the group under which the d'Alembertian operator in eq. (jl]) is invariant when 
c 2 > 0. Let s 2 = x 2 + y 2 — (cz) 2 and let r\ 6 (— oo, oo) be a real- valued hyperbolic angle, and 
let (p G (— 7r, 7r] be the usual azimuthal angle. Then the following parametrizations cover the 
"future," "past" and "absolute elsewhere" of the origin, respectively: 
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s sin ( 
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Suitable homogeneous polynomials that satisfy eq. (j3J) in these three regions and belong to 
irreducible representations of £0(2, 1) are 

(s 2 f 2 i™Y itm (i V ,<t>) (s 2 y/ 2 i m Y e , m (*-ivA), and (s 2 f 2 i m Y e>m ^ - i V , 0). (18) 

The tidal response of a coreless incompressible planet is made up of linear superpositions 
of these: for example, eq. H16c[) is proportional to s 2 Y 2 fi(iri, 4>) i n the "future." In order 
to accommodate the boundary condition vr = at the surface of a finite core, one would 
like to add to the polynomial solutions (1161) a suitable linear combination of homogenous 
solutions of negative degree. The negative-degree solutions would be larger near the core 
than near the surface, and therefore could "patch up" the boundary condition at the core 
without much spoiling the boundary condition at the surface. The obvious negative- degree 
counterpart to the first of the functions eq. f|T8|) . by analogy with the case of Laplace's 
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equation, is (— s 2 )~^ +1 ^ 2 i m Ye >m (ir], <p). While this is indeed a solution of the wave equation, 
it is unfortunately singular where the core intersects the "light cone" s 2 = 0, i.e. at the 
critical colatitudes 6 = (3 and 6 = n — (3 [eq. ffTUl) ]. So it seems that we cannot superpose a 
finite number of finite-degree homogeneous solutions to construct a smooth response. 



3.2. WKB scattering calculation 

3.2.1. Nonspecular reflection from a planar boundary 

To see how short inertial waves may be produced from long ones, consider the reflection 
of an incident plane wave, 

Mx,t) = A^e*-"-™ , 

from a planar wall with normal n pointing away from the fluid. The reflection is generally 
not specular, and the incident and scattered wavelengths differ, unless n is parallel or per- 
pendicular to (I. These are consequences of the anisotropy of the dispersion relation Q and 
would also hold, with some quantitative changes, were buoyancy important. 

If the wall is fixed in the unperturbed frame of the fluid and sufficiently rigid so that 
any transmitted wave is negligible, then energy conservation requires a reflected wave. This 
is enforced by the boundary condition, which we have taken to be n • v = 0, but any other 
non-absorbing boundary condition would lead to the same relationship between the incident 
and reflected wavevectors as the one we are about to derive. 

The scattered (outgoing) wavevector k' is determined by two conditions. First, in order 
that the incident and scattered wave have the same relative phase at all points along the 
wall, as required by the boundary condition that connects them, it is necessary that k and 
k' have the same components parallel to the wall: that is, n x (k — k) = 0, or equivalently 

k! = k'j_n + (k- nfi-k) = k'j_h + k\\ , (19) 

where k' ± remains to be determined. Similarly, in order that the relative phase be constant 
in time, the two waves must have the same frequency, uo. Substituting from eq. (1191) for k' 
into (k'foo 2 = (2fi • k') 2 leads to 

[(2fi • n) 2 - uj 2 ] {k'J 2 + 8(0 • fe||)(n • h)k' ± + (2(1 ■ fey) 2 - ui 2 k\ = 0. (20) 

It is clear that one of the roots of this quadratic equation must be the known solution 
k' ± = k± representing the incident wave, and therefore the product of the distinct roots is 

u 2 kl - (2(1 ■ k\\) 2 
k - k '-= J-(2»-n)" ■ < 21 > 
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The denominator in eq. (|2ip may vanish when the numerator does not: the reflected wavevec- 
tor then becomes infinite and normal to the boundary. We believe that the change of wave- 
length upon reflection, which is a direct consequence of the inertial-wave dispersion relation, 
underlies the singular behavior observed by OL04 in their numerical calculations. 



3.2.2. Scattering from a spherical core 

Now we apply this idea to scattering of the "equilibrium tide" from the spherical core. 
As noted in §2.31 the equilibrium tide is radially nodeless, and in this sense nonwavelike. 
However, the functions fTTOl) do have nodes in angular directions: the wavenumber parallel 
to the surface of the core of these quadratic functions is roughly 

2 m 

fey « ±— e e ± — — — -e^ = kgig + k^ (22) 
rC c ri c sin o 

at colatitude 9. Since |fei|i2 c | > 1> we may expect that WKB should be applicable to an 
outgoing wave whose radial wavenumber k' ± is real and ^> 1/R C . We cannot make use of 
eq. (}2"Tj) because the incident component k± is ill-defined and the equilibrium tide does not 
satisfy the WKB dispersion relation. But equation (1201) should still be applicable to the 
short-wavelength outgoing wave. It can be seen that the coefficient of k' ± vanishes at the 
critical latitudes 9 = 3 and 9 = n — 8, where 8 is defined by eq. ( TIP}) . 



The large root of eq. fl2"01 can be found by balancing the terms in k' ± and (k' ± 



\2. 



k' 



8(12- fey ){n-h) k e sin29 



R ~ ,,2 ^OO.^A2 



UJ 



(2fi-n) 2 cos 2 9- cos 2 /? 



k e csc(8 — 9) near 9 = 8 

v ; (23) 
ke csc(/3 + 9) near 9 = tt — 3. 

Let us now determine when the large root represents an outgoing wave. Since \k'^\ 3> fen, we 
may write k w sign(fc^)e^ + |A;^| _1 fe|| through first order in k\\/k' R , whence from eq. (11 II) . 



F gr « 2fisin/3sign(cos^) {jjfyZR - Y* e ) ^ 

near the critical latitudes. The radial component of this is positive when sign(/c0) = 
sign(cos#), so this becomes the condition for the outgoing wave: kg > at the northern 
critical latitude, and kg < at the southern one. Also, it can be seen from eqs. (1241) & 
that the group velocity is predominantly latitudinal and directed away from the critical 
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Fig. 1. — Scattering of a long wave into a short one at the core near the northern critical 
latitude. Allowed directions for k make angles (3 or ir — (3 with ft, as shown by "light cone." 
Large diagonal solid arrows (red in online version) are outgoing fc out . These become infinite 
and change sign across critical latitude (thin dashed line). Smaller solid arrows (blue online) 
show group velocity Vg )OU t, whose magnitudes are oc k~ x . Heavy dashed arrows show the 
incoming fc in and its associated group velocity. These are not to scale: k lTL <C k ont and 
Vg,m 3> ^g.out- In fact, all of the red arrows should have latitudinal components that are 
nearly the same in magnitude and sign, w kg in > 0. 
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latitude: that is, towards the pole on the poleward side, and toward the equator on the 
equatorial side. 

Next, we determine the amplitude of the outgoing wave by matching it to the positive- kg 
component of the incident wave via the boundary condition Br • (v in + v out ) = 0. The radial 
velocity Vr is related to if) by eq. (fT5|) in cartesians, and in spherical polars [cf. eq. (j3J)] 

1 



9 „ 9 „. dib sin 26 dib m cos 3 , 



(25) 



iuj sin 8 

For the outgoing wave, this reduces at the critical latitudes to 

| kg | R c sin — im 
2\ IR C sin p 

The contribution from dRib is finite and oc kg despite the fact that k' R diverges, because of 
the factor (cos 2 # — cos 2 /?) in front of dib out /dR. For the "incoming" equilibrium tide, the 
expression is similar, except that dnif> in doesn't contribute because dip in /dR is finite at the 
critical latitude: 

— | kg | R c sin f3 — im 
2\ IR C sin p 

Combining the last two equations, we have the following relation between the long- wavelength 
"ingoing" tide and the short-wavelength outgoing wave: 

A^JM^tl^A^ (26) 
\kg\R c sin p + im 

Thus, the boundary condition £i?,out = — £,R,m at the core leads to |A out | = \A in \. 

The relationship between the amplitude A- m of the "incident" component of the equilib- 
rium tide and the tidal potential $i is given by equations (flEl) . For definiteness, consider a 
synchronization tide exerted on a planet in a circular orbit, so that 

2n 

m = 2sign(n — Q), and u = 2\n — Q\, where n = — — (27) 

-'orb 

is the mean motion. In this case, the relevant component of $i becomes 

$1 _ _M^Ml^ S in 2 g e ^-- * . (28) 

The factor / represents the ratio $i/3>i, e xt of the total perturbing potential to that part 
which is exerted by the companion. In general, / ~ 1 + 2k p at the surface, R = R p , where k p 
is the apsidal- motion constant, equal to one half the planetary Love number. For a constant- 
density coreless planet, / = 5/2. Expanding sin 2 # and extracting the coefficient of exp(2?#) 
leads to 

9/ u(u + mtt) GM*Rl r , 

A ™ = -T^^1TT; " H for m =2ande = 0. 29 

16 87rGp a A 
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3.2.3. Energy flux and power of the scattered waves 

Since |VWt| = |An| equations ffT2l) . (jUj) , and (1291) imply that the radial component of 
the energy flux in the outgoing short waves at the core is 

~ ~ V 32M p a 3 J V^^ P c { ' 

for the \m\ = 2 synchronization tide, i.e. for the tidal response of a planet in a circular orbit 
with aligned but nonsynchronous rotation. Expressions for u and m in terms of the orbital 
and rotational periods were given above in §3.2.21 As it stands, eq. fl30|) is valid only for 
a const ant- density planet (so that / = 5/2), and only for the waves launched sufficiently 
close to the critical latitude so that we may take \k \ ~ \k' R \ 3> kg, k§. We have taken 
\k$\ ~ 2/R c and k^ = m/R, as appropriate for the response ( jj~6l) to an \m\ = £ = 2 tidal 
potential. Notice that the large radial wavenumber has canceled between the wave energy 
density (j 12 aft and radial group velocity ([21]), so that the flux is approximately independent 
of latitude, provided that 9 is sufficiently close to /3 — cos~ 1 (uj/2Q) or to 7r — (3. Let us 
therefore integrate over latitudinal bands of width 2A6 centered on both critical latitudes, 
with A6 <C [3, tt — (3, to obtain the total mechanical power carried radially outward by the 
short waves launched within these bands: 

E(A0) w 8nR 2 c sm{3A6 x T R 

« 2.(^) (1)(u+mn) ^_^j P R*A9. (31) 

Although eq. (13T|) has been derived for an incompressible, const ant- density planet, it 
can be generalized to an isentropic compressible body with standard approximations for 
the equilibrium tide. In our approach where the core is regarded as a perturbation to an 
otherwise homogeneous body, we continue to calculate the equilibrium tide as if the core 
were absent. What is needed for the scattering calculation is the radial displacement of 
the equilibrium tide at the coreQ ^(R c ) For a const ant- density body, the radial strain 
a = £ R (R)/R is independent of radius R, so 

t R {Rc) ~ T P U {Rp) ~ ~R~ P ~gV~ ~~ f {~) c ' ( } 

where a is the semimajor axis of the orbit, g p = GM P /Rp, and the surface displacment has 
been evaluated from (El) with neglect of ip on the righthand side as before, i.e. (,r(R p ) ~ 



1 It is assumed here that £r and <!>! are proportional to the same function of (t, 9, <f)) as the perturbing 
quadrupolar potential $i iG xt of the star, so this dependence will be taken as read 
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—Q>i(Rp)/g p . We assume that the relation (132]) between the surface displacement of the 
equilibrium tide and the disturbing potential holds in general, provided that / is interpreted 
as 1 + 2k p for the relevant value of the apsidal motion constant k p . The latter is k p « 0.260 
for an n = 1 Emden polytrope, which roughly approximates a Jovian planet, as compared 
to k p = 3/4 for a constant-density body (n = polytrope). However, since the radial strain 
of the tide is not constant with radius in general, we must multiply / in eq. (1321) by a factor 
cr(R c )/a(R p ) to obtain the correct displacement at the core. If the core is sufficiently small, 
R c < R p /2, then a(R c ) « a(0). 

It is not immediately clear how to estimate cr(0) easily. One possibility is to assume that 
which represents the radial displacement of the fluid by the equilibrium tide, is the same 
as the radial displacement of the equipotential surfaces between the undistorted and distorted 
states; then it can be shown by integration of the Radau equation from which k p is obtained 
JSchwarzschildl Il958h that a(0)/a(R p ) « (1 + 2A; P )' 1 « 0.658 for the N = 1 polytrope, 
so that the factor in square brackets in (1321) reduces to unity for N = 1. Alternatively, 
if one evaluates from the summed zero-frequency response of all the normal modes, 
taking into account their overlap with the quadrupolar perturbing potential $i iex t, then it 
can be shown that a(0)/a(R p ) ~ 0.813 at N = 1. By far the largest response is that of 
the I = 2 fundamental mode, for which £r is nearly linear in R. These two estimates of 
a(0)/a(R p ) differ because, even for an originally spherical and nonrotating body, the radial 
displacements of the fluid and of the equipotentials need not coincide except at the surface. 
There are many possible displacement fields that can be compatible with a given distortion 
of the density and potential fields, and it is not possible to choose among them on the basis 
of the continuity equation alone without some auxiliary constraint. It can be shown that 
in a stratified region, the appropriate constraint is V • £ = even for a compressible body: 
the density and pressure of fluid elements is not disturbed by the equilibrium tide in the 
stratified case. OL04 adopted this constraint. For an isentropic region in a nonrotating star, 
the appropriate constraint is V X $, = rather than V • $, = since vorticity is conserved; 
in a compressible body, this leads to a di fferent pattern of fluid disp l acements for the sam e 
distortion of the density and potential (jGoodman fc Dickson 1998 ; Terquem et al. 1998 ). 
Thus our second estimate (the one yielding 0.813) was calculated under the assumption that 
£ = V% for some scalar function \. Since our bodies rotate, however, the axial component 
of vorticity does not vanish and therefore V X £ = may not be the correct constraint @ For 
the purpose of estimating the order of magnitude of the tidal Q that results from scattering 



2 In the strict low- frequency limit where uj 2 is small compared with 40 2 , and not just small compared to 
w dyn' the vorticity would remain axial under the equilibrium tide by the Taylor-Proudman theorem, so that 
the poloidal components of £ would be curl-free. 
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by the core, the differences among 0.658, 0.813, and unity (the correct result for N = 0) are 
not important, and we take this as an indication that the true value of a(0)/a(R p ) is also 
sufficiently close to unity. 

The frequency- dependent factor u(u + mfl) in eq. (TSTI) becomes 8\n — Q\n 2 with use of 
(|27|) . Finally, the relevant density is the density of the fluid at the surface of the core, p(R c ). 
With these modifications, the wave power (13T]) generalizes to 

©'(t)'**^- (33) 

The tidal torque, or more precisely, the rate of increase of the angular momentum carried 
by the waves, is T — E x m/u — E / (n — fl) . Thus according to eq. (1331) . the torque has the 
same sign as n — Q — meaning that subsynchronous planets spin up and supersynchronous 
ones spin down — but is independent of the magnitude of the departure from synchronous 
rotation. This is probably not true if the departure is very small, \n — 0| <C n: in that limit, 
the critical latitudes converge upon the equator from both sides, whereas the approximations 
used to derive eq. (|33|) implicitly assume that these latitudes are well separated from one 
another. 



E(A6) 



i oi 2 97r 
n — si n — 
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3.2.4- Dissipation of the short waves, and the tidal Q 

Waves that do not dissipate appreciably before returning to the region from which they 
are launched (perhaps after multiple reflections between the outer boundary and the core, 
with changes of wavelength at each reflection) must be treated as global normal modes. 
Secular input of energy and angular momentum to nondissipative global modes would occur 
only at exact resonance with the tide, which would almost never occur for modes of finite 
wavelength. Furthermore, the angular momentum carried by the waves is not transferred 
to the mean flow, and therefore does not alter the angular velocity Q, until those waves 
dissipate. 

Thus, in order to estimate Q, we must consider the dissipation of the waves. Unfor- 
tunately, this is a complicated issue. More than one process may be important, depending 
upon aspects of the planetary structure and transport processes that can justifiably be ne- 
glected in calculating the wave excitation. Nevertheless, some dissipative processes can be 
ruled out, and a rough upper limit on Q as a function of the core radius can be obtained. 

First of all, the narrower the width A8 of the latitudinal bands around critical latitudes 
that we consider, the shorter is the wavelength of the waves launched within those bands: 
eq. fl23l shows that A ~ ttR c A6. All of the obvious dissipation mechanisms become more 



co- 



efficient as wavelength decreases. The rate of viscous dissipation, for example, scales as vk 2 , 
and therefore oc |A#| -2 . Furthermore, the propagation time from the core to the surface of 
the planet — or rather, to the upper boundary of the convection zone, since in a realistic hot 
Jupiter there must be a stratified region near the surface — scales oc A oc A9. [Eq. (T23J) says 
that the radial component of the group velocity starts out at V gVy R ~ QX 2 /R C . But this is 
because the rays emanating from near the critical latitude are almost tangential to the core. 
The rays follow straight lines in the meridional plane, and A is approximately constant along 
them, so that a more representative value of the group velocity for the purpose of estimating 
the propagation time is |V gr | ~ Q\ ~ QR C A9.} Hence the number of viscous dissipation 
times per transit time scales oc |A#|~ 3 . 

Nevertheless viscous dissipatio n is probably negligible , as c an be seen by very rough 



order of magnitude considerations. iGoldreich fc Nicholson! (119771 ) estimated that turbulent 
convective viscosity acting on the equilibrium tide, which has a "wavelength" ~ R p , would 
yield Q ~ 10 13 because of suppression of turbulent viscosity by a factor (ut c )~ 2 , where r c is 
the turnover time of the largest eddies. (When this argument applies, the Q based on the 
laminar viscosity should be even larger.) In our case, the suppression factor would not be 
quite so small because the tidal period is a few days even for a substantially nonsynchronous 
rotation, rather than 5 hours as for the Jupiter-Io system. So by Goldreich and Nicholson's 
reasoning, Q ~ 10 11 for the equilibrium tide in our case. The short waves have the same 
period as the tide itself, so the suppression factor is the same for them, but their damping 
rate is increased by a factor ~ (kR) 2 ~ IA6 1 ! -2 , and allowing for their transit time between 
the core and the surface, we may conclude that the Q value due to turbulent convective 
damping of short waves should be at least 10 n |A#| 3 in hot Jupiters. We have stated this 
as an inequality because the core is small and therefore somewhat inefficient at scattering 
the equilibrium tide; this will be made more quantitative below, but for now, note simply 
that the wave power (T3HT) is oc R b c . To be astrophysically relevant, Q should be of order 10 5 
to 10 6 . Therefore only waves within A9 < 10~ 2 radians of the critical latitude could damp 
effectively by this mechanism. But as we will soon show, other mechanisms exist that can 
damp the waves launched farther from the critical latitude, and since the wave power fl33l) 
is proportional to A9, these mechanisms give a smaller Q. 

An important source of dissipation for short inertial waves in hot Jupiters is escape from 
the convection zone, where most of the mass of the planet resides, into the stratified radiative 
zone near the surface, whose existence is guaranteed by strong illumination from the host 
star. In the radiative zone the waves convert into g modes (more properly, Hough modes) 
that are supported primarily by buoyancy rather than Coriolis forces. They are then subject 
to damping by radiative diffusion because they perturb the temperature and entropy profiles. 
Radiative diffusion damps more efficiently than viscosity because of the much longer mean 
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free path of photons compared to molecules or ions, and it is all the more efficient because 
the wavelength in the radiative zone shortens by a further factor ~ Q/N ~ 10 -2 , where 
N is the Brunt- Vasala frequency there. The upshot is that an outgoing short inertial wave 
that penetrates the radiative zone will almost certainly damp before returning to the core. 
Details will be given by Eric Johnson in a forthcoming paper. Here we simply note that 
penetration is not possible between the poles and the critical latitudes, because the Hough 
modes are evanescent there. Using the fact that the group velocities lie at angles (vr/2) ± (3 
with respect to the polar axis in the meridional plane, one can show with a little trigonometry 
that penetration into the radiative zone can occur at the first encounter only if 



R c 



cos A6 < 



1 (poleward rays) 
(equatorward rays), 



(34) 



where R s ~ R p > R c is the radius at the convective-radiative interface. A poleward ray is one 
that is launched between the critical latitude and the pole, so that it starts out toward the 
rotation axis, etc. Neither inequality in ( 1341) can be satisfied when u < i.e. (3 > 7r/4. 

In such cases the waves are fully reflected at their first encounter with the radiative zone and 
return toward the core. But the ray may enter the radiative zone on a subsequent encounter. 

Paradoxically, the easiest damping mechanisms to predict with confidence by analytic 
means may be the nonlinear ones. It is reasonable to assume that any inertial wave whose 
velocity amplitude satisfies v > u/k will quickly damp by some combination of Kelvin- 
Helmholtz instabilities or three-mode coupling to even shorter-wavelength daughter modes, 
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lengthscales comparable to the wavelength and timescales comparable to the wave period. 
The dimensionless nonlinearity parameter kv/u diverges rapidly toward the critical latitude. 
From eqs. (fl2"j) . \v | 2 = 4jF H / pV gt R] substituting then from eqs. (1241) and (130|) . and generalizing 
the latter in the same way that we turned eq. (j3~il into eq. (J33l) . we find that 



kv 



2V8 



1 + 2fc c 



<j(R c ) lo + mVt 
cr(Rp) \ uj — mVt 



M * R * (MY 



(35) 



The contents of the square brackets above are close to unity. Therefore, nonlinear dissipation 
will dominate within latitudinal bands of halfwidth 



A6 



ni 



M*Rl 
M p a 3 



1/2 
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radians. 



(36) 
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Though small compared to unity, this is indeed larger than the generous estimate made 
above for the latitudinal distance within which turbulent convective viscosity might be im- 
portant. We regard this as a lower bound on the value of AO within which the waves are 
able to dissipate, since other mechanisms — especially escape into the radiative zone — may 
contribute. 

The tidal Q = uAE/E, where AE is the maximum potential energy associated with 
the time-variable part of the tidal distortion. The instantaneous gravitational energy of the 
equilibrium response to an applied quadrupole $i jext (-R, 9, 0) = KR 2 P 2 (cos'y) is k p K 2 R^/G; 
here K is a constant and 7 is the angle between the symmetry axis of the potential and the 
point R. For the present case of a synchronization tide, K = — y^3/2GM^/a 3 [eq. (IA.1I) ]. so 



3 7 GM 2 Rl 

2 V 



ae -/,-,. f ~ " 



(37) 



Using eq. (T2TI) & (1531 and taking k p « 0.26 and p(R c ) ~ (7rM p /AR 3 ) as appropriate for the 
apsidal motion constant and central density of a homogeneous N = 1 polytrope, we have 
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A6> 



0.1 rad 



-1 



(38) 



In the final line, we have taken the square brackets equal to unity. This result is some two 
orders of magnitude larger than is typically assumed for hot Jupiters but is obviously very 
sensitive to the assumed core radius. The fiducial value R c ~ 0.2i?j is a crude estimate based 
on the assumptions that M c w 0.1M P and that the core must be roughly twice as dense as 
its immediate surroundings at the same pressure because it has roughly twice the molecular 
weight per electron. Here P is the period of the orbit (P = 2k /n), not the period of the tide 
(-Ptidc — 7r/\n — Q\): Q turns out to be independent of the latter. The Q for circularization of 
a slightly eccentric but synchronous orbit is nearly the same for this mechanism if the same 
AO applies (see the Appendix). 

Values for Q quoted in the literature are often n ormalized by th e ti dal energy of a 



homoge neous body of the same mass and radius; see iMardlingi (120071 ) & iFabrycky et al. 
(l2007bl ) for discussion of this point. Since k p = 3/4 rather than w 0.26 for a homogeneous 



body, the values (1381) and (IA7j) should be roughly tripled to conform with that convention. 
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4. Discussion 



Up to this point, we have had little to say about the work of Yanqui n Wu on the 
excitation of inertial waves in coreless but compressible isentropic fluid bodies (Wu 2005al lbh . 
Here we explain why we believe that her calculations, though technically superb, are not 
applicab le to tides in real pla n ets or stars. Some of these criticisms may apply also to the 
work of llvanov fc Papaloizoul (120071 ) , who adopted a "low- frequency" approximation that 
appears to be physically equivalent to Wu's. But we focus on Wu's work because of the 
admirable clarity of her exposition, and because she is concerned with resonant ex citation 
rather than the quasi- impulsive excitation analyzed by llvanov fc Papaloizoul (120071 ) . 



We then go on to discuss whether the rock-and-ice core of a jovian planet can be regarded 
as rigid, or even solid, and the consequences for the production of inertial waves if it is not. 



4.1. Previous tidal calculations for coreless isentropic bodies 

The salient claims by Wu that we address are the following: 



1. The overlap integrals between the short- wavelength inertial modes and the perturbing 
tidal potential diminish as a negative power of the number of radial or latitudinal 
nodes, rather than exponentially, even for unperturbed radial density profiles that are 
smooth apart from a power-law convergence to zero at the surface. 

2. For many radial nodes, the overlap is concentrated toward the surface; this is where 
the excitation mainly occurs. 

3. The overlap integrals vanish for a constant-density, incompressible body. 



Wu suggested that the tidal coupling could be enhanced by a density discontinuity associ- 
ated with a first-order phase transition within the convection zone, perhaps at the interface 
between molecular and metallic hydrogen. This is probably true in principle, but we ad- 
dress here only the claims made for strictly isentropic bodies, thereby excluding first-order 
phase transitions because of the entropy jump associated with latent heat. There does not 
appear to be a consensus as to the first-order nature o f the molecular-to-metallic transition 
in hydrogen (IMilitzer et al.l 120081 ; ISumi fc Sekinol 120081 . and references therein) . 



We agree in part with the third of the ennumerated claims above but take issue with 
the first two. Our counter-arguments are mainly these: 
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1. By selectively neglecting compressibility in some places while retaining it in others, and 
especially by oversimplifying the low-frequency limit of the outer boundary condition, 
Wu has created a singularity at the boundary, though the dynamics there should 
actually be smooth if the boundary is free. We suspect that the tidal coupling she 
calculates is due mainly or entirely to this singularity. 

2. For barytropic bodies [P = P{p)], albeit with an inconsistent treatment of their grav- 
itational potentials, there exist tidal responses that completely lack short-wavelength 
components and that can be exhibited in closed form. The fluid displacement and 
enthalpy are independent of adiabatic index, and so are the same for a compressible 
as for an incompressible body. 



We now expand upon these last two points. 

After deriving the equivalent of equation (TjJ, Wu argues that the lefthand side, which 
contains the sound speed c s in the denominator, can be neglected on the grounds that 
compression of the fluid is very slight for modes that are both short in wavelength (kR ^> 1) 
and low in frequency (uj 2 <C c 2 k 2 ). This term comes directly from the time derivative in 
the continuity equation, so dropping it is equivalent to replacing the continuity equation by 
V-(pvi) = 0. But accord ing to her formalism, the tidal forcing is explicitly shown 



by eq. (7) of IWul (j2005bl ). so the tidal coupling itself should vanish in the limit c~ 2 — > 0. 
Indeed, her Appendix B makes this explicit. However the neglected term is singular at the 
free surface, because c 2 — > there. In the full equation (j3J, the singularity involving c~ 2 
is balanced by another singularity involving oc p _1 Vp, where p is the unperturbed density, 
Using our equations (jSj) & (jlj), it can be seen that the condition under which the two 
singularities cancel one another is precisely the free boundary condition (jSJ), whose physical 
interpretation is the vanishing of the lagrangian enthalpy perturbation. 



If, following (IWul l2005al . §2.1 & §4.1), one neglects the righthand side of eq. (JH) but 
retains the gradient of the unperturbed density profile on the righthand side, then the only 
possible well-behaved modes of the system are those for which the normal components of 
the fluid velocity and displacement [£ = v/(—iu)] vanish at the boundary. Evidently, Wu 
believed that taking n ■ £ = at the boundary is an acceptable approximation because (i) 
the inertial modes don't move the boundary very far, and (ii) the density vanishes at the 
boundary anyway. The modification to the boundary condition affects the structure of the 
mode not just at the surface, however, but down to depths comparable to that of the first 
radial node of ip below the surface. Let h be half the nodal depth, so that ip(h) ~ t/>(0), 
where ip(Q) is the surface value of the eigenfunction. Denote the unperturbed density at 
this depth by p(h), and the sound speed by c s (h) = a/ gh(N + l)/iV, assuming a polytropic 
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equation of state P oc p( N+1 )/ N . The eulerian density perturbation at this depth is 

Pl(/l)=p(/l) ^)"ivTT^- (39) 

(With Wu's definition of ip, there would be a factor of u 2 here.) On the other hand, since 
p(R) oc (Rp — R) N near the surface, the contribution to p\ = — V- (p£) from the surface 
displacement, which Wu neglects (as do Ivanov & Papaloizou), is 

-Z ■ V P\R=R P -h ~ — ^— &(0) ~ — ^ — , (40) 

where we have evaluated £r(0) from the free boundary condition (j3J) in the last step, taking 
$1 = as appropriate for a free mode of oscillation with negligible self-gravity. Evidently, 
the neglected contribution (HO"]) is comparable to the total fl22D, an d therefore is not negligible 
within the first node. 

In Wu's powerlaw-sphere model, where 

p oc (R 2 max - R 2 f (41) 

the nodal depth is oc R/n, except near the critical colatitudes 6 = cos _1 (±c<j/2f2) (the 
"singularity belt" in Wu's parlance) where it scales oc R/n 2 . Here n = n% + ri2 in terms of 
Wu's modal indices rij > 0, which are roughly proportional to the WKB wavenumbers of the 
inertial modes: that is, n ^> 1 for shortwavelength modes. Wu finds that the tidal forcing, 
which is proportional to the overlap integral between the perturbing potential $i and the 
modal eigenfunction, scales with this index approximately as n~ 2l3 ~ l . This is consistent with 
the idea that the excitation occurs within the first node from the surface (and probably also 
within the singularity belt) since the fraction of the planetary mass within depth 2h of the 
surface for the density profile flUD is AM/M oc (2/i) /9+1 . 

It is true that the horizontal components of the velocity and displacement near the 
boundary are larger than their radial components by a factor ~ ui\ yn /Vt 2 ^> 1, where c^dyn = 
(3g/i?) 1//2 is the dynamical frequency of the planet. So it is likely that the errors in modal 
energies and eigenfrequencies caused by the approximate boundary condition ^ R = are 
slight. But the overlap integrals are of a higher order of smallness in (Q/udyn) 2 , so that their 
relative errors could be large. 

The discussion so far does not make clear the sign of the error (if there is one) : perhaps 
the tidal coupling would be large r with the exact bound ary condition. The following model 



system, which is borrowed from iGoodman et al.l (119871 ). suggests that the error is in the 



direction of overestimating the coupling. 
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Let the equation of state again be polytropic, with p oc w N , P oc w N+1 , and c 2 = w/N, 
where w = (N + l)P/p is the enthalpy, which remains finite in the incompressible limit 
iV — > + . To match (|41|) . the unperturbed enthalpy should be 

w(x, y, z) = w(0) (l - x2 + y 2 + z2 \ f (42) 
and N = (3. Hydrostatic equilibrium in the corotating frame (where v — 0) requires 







w + $ - ifi 2 (x 2 + y 2 ) 



(43) 



where $ is the unperturbed potential. Since w and the centrifugal term are quadratic 
functions, $ must also be such a function. In this regard, the model differs from the one 
considered by Wu. In effect, she calculates the potential due to the density profile (j4"Il) 
from Poisson's equation and uses this to determine the enthalpy. Consequently her pressure 
profile is not simply a power of the density profile, and I\ = (d In P/d In p)g is not constant 
in her models, at least not for (3 > 0. However, she argues at several points^] that the tidal 
forcings she calculates depend, at least for a smooth density profile, only on the behavior 
near the boundary, where it is indeed approximately true that P oc p(^ +1 )/^. As a matter 
of fact, because of her neglect of the c~ 2 term in eq. (JH), the equation of state doesn't enter 
her calculations of the normal modes: only the density profile does, which could result from 
many isentropic equations of state paired with an appropriate background potential. The 
overlap integrals do depend upon the equation of state via the sound speed, but as long as 
c 2 approaches zero linearly near the boundary, it is hard to see how the overlap integrals for 
short- wavelength modes could be sensitive to the full functional forms of c 2 (i?), and therefore 
of $o(-R) ; if they are excited near the boundary. For these reasons, it does not seem crucial 
that the unperturbed potential be fully consistent with the mass distribution. 

After elimination of the density in favor of the enthalpy, the linearized equations become 

-iuv + 2f2 x v + Vwi = -V$i, (44a) 
-iuwi + v • Vw + N^w V- v = 0. (44b) 

Now suppose that the tidal potential is quadrupolar: that is, a homogeneous and harmonic 
quadratic polynomial in (x,y,z); for definiteness, 

$ 1 = ^(x + ^) 2 e-^ (45) 



For example, in after eq. (C7)): 'The results only depend on the boundary behavior of 



/(0) as long as it is sufficiently smooth. This explains why models with different polytrope representations 
(p oc [1 — (r/Ryf or p oc p 1+1 ^) give rise to essentially the same overlap integrals." 
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where A is a constant, and the exponential factor will be taken as read hereafter. Then since 
w is also a second-degree polynomial, eqs. fj44l) can be satisfied by taking the components of 
v to be polynomials of the first degree, and vj\ of the second degree. After some algebra, 

v x = -iv y = a(x + iy), v z = 0, w x = —(x + iy) 2 , 

, , R 4w(0) 

where a = — ; — r . . . ^ — , and B = — — — a . 4b J 

Since V-t> = 0, the equation of state (i.e. N) doesn't enter the solution (1461) . Also, the 
relative vorticity Vxt> = 0, so this solution has the same total vorticity 2f2 as it would have 
in the absence of the tide and therefore might be the solution of an initial value problem in 
which the tide was "turned on" slowly. The denominator in the expression for a vanishes at 
to = — ± (lUq — fi 2 ) 1 / 2 , where ujq = 2^Jw{ti)/R max is the dynamical frequency of the model 
and therefore presumably 3> VL. At the unperturbed surface where w — 0, eq. (I44bl) reduces 
to the free boundary condition that the lagrangian enthalpy perturbation vanishes. 

These details aside, the important point is that the tidal response is entirely long- 
wavelength for any polytropic index when a free rather than rigid outer boundary condition 
is used, at least in this idealized coreless model, which uses a nonselfconsistent but smooth 
unperturbed potential, and at least for a nonsynchronous body in a circular orbit. Short- 
wavelength inertial modes are not tidally forced even though the fluid is compressible. 



4.2. Rigidity of the core 

OL04 assumed the core to be solid, and therefore impenetrable by low-frequency but 
short-wavelength disturbances such as inertial waves, but sufficiently plastic as to comply 
with the large-scale equilibrium tide as if it were fluid. Here we re-examine the strength 
and solidity of the core. In agreement with OL04, we find that the elastic strength of even 
a solid core would be negligible as regards the equilibrium tide. However, we evaluate the 
equilibrium tide differently than OL04, and we allow for the density contrast between the core 
and its immediate surroundings. Furthermore, we estimate that the core is most probably 
fluid rather than solid. 

It is presumed that the cores of Jovian planets consist of elements heavier than hydro- 
gen and helium, more specifically of some combination "rock" (refractory m inerals such a s 



silicates and iron) and "ice" (molecular species such as H 2 0, CH 4 , and NH 3 ) (iGuillotl 120051 ). 



To support shear stress, these materials would have to be in a solid phase. The pressure at 
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the surface of the core is comparable to the central pressure of an N = 1 polytrope: 



- 4xll UJ UJ dyncm • (47) 

or ~ 40Mbar. The temperature is more difficult to predict. Jovian planets are supported 
mainly by degeneracy pressure, so the central temperature has only a modest effect on the 
planetary radius. The temperature depends upon the planet's age and rate of cooling, which 
in turn depend upon on uncertain opacities in the envelope, not to mention the possibility 
of internal heating by tides. Present estimates are in the range 

T c ~ 2-4 x 10 4 K (48) 

based on stand ard models fit to radii of transiting planets and the estimated ages of their 



host stars (e.g. lArras Bildstenl 120061 ) . 



The pressure (I47j) is large compared to bulk moduli of common refractory materials at 
room temperature — e.g. K = 1.7 and 1.0 Mbar for iron and silicon, respectively — so rocky 
cores will be compressed to de nsities > 20gcm~ 3 [ we de rive this number from a generic 



equation of state for "rock" by iHubbard fc Marleyl fjl989f )] . i.e. a factor 3-10 times larger 



than their densities at atmospheric pressure. Under standard conditions, the shear moduli 
of such materials are comparable to their bulk moduli (pp e ~ Psi ~ 0.8 Mbar). Under 
compression, the bulk moduli rise more quickly than the shear moduli because the former 
is associated with the degeneracy pressure of the electrons, whereas the latter is a Coulomb 
effect having to do with the ion lattice; for very large compression factors, one therefore 
expects n/K oc (p/po) 1 / 3 . With these scalings, we can compare the elastic energy of the 
core under its distortion by an equilibrium tide with the corresponding gravitational energy. 
Approximating the core by a constant density p c , we have 



AF K GM* 2 3pM 2 

A£ elas _ 25pfl c 
A£ gTav ~ 4GM c p 



c 



(49) 



where a is the radial strain £ r ,max/-R c caused by the quadrupolar tide. The elastic modulus in 
the core's compressed state should be p ~ (p/p )^ 3 po ~ 2p , where the subscript "0" refers 
to atmospheric conditions, hence no more than a few Mbar, whereas GM c p c / R c m 200 Mbar 
if M c = 30 M e and p c = 20gcm~ 3 , which is of course comparable to the pressure (1471) . So 
the elastic energy of the equilibrium tide in the core is only 5% of the gravitational energy, 
and therefore the large-scale tidal distortion of the core should be nearly the same as for a 
fluid. However, elastic strength could still be enough to prevent the propagation of inertial 
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waves inside the core, because the tidal period is long compared to the crossing time of an 
elastic wave in the core, which is of order half an hour for the numbers above. 

The discussion so far has been based on a solid core, but a liquid one may be more 
likely. Melting of an ionic lattice tends to occur at T = (Z eff e) 2 /Yj£;T ~ 10 2 , a dimensionless 
measure of the relative importance of Coulomb to thermal ener gies; here Tj ~ (3/3/47rm ? ;) 1 / 3 



is the mean distance between ions in terms of their mass, rrii (e.g. IShapiro fc Teukolskylll983 
and references therein). The question is what to use for the effective charge Z e $ governing 
the ionic interactions. For silicon at its reference density of 2.33 gem -3 , for example, this 
formula would predict a melting temperature ~ 2 x 10 5 (10 -2 r)~ 1 K if one were to take 
Z cfi = 14, the full charge on the nucleus. In fact, most of that charge is shielded by electrons 
whose orbits are much smaller than r^, so that Z e ^ w 1 is a more reasonable choice; this 
leads to QQQT^ 1 K, which is comparable to the actual value, T m = 1693 K. As far as we 
know, there is no experimental measurement of the melt temperature near 44Mbar, but if 
one simply scales it from 1 bar by the the reciprocal of the inter-ionic distance, assuming 
that the Coulomb interaction is characterized by a constant Z e g throughout this range, then 
T m @20gcm -3 ~ 3500 K. Despite the crudeness of this argument, it therefore seems likely 
that "rock" should be molten at the much higher temperatures in eq. (j480 . 

To recap, the core is probably fluid, and even if it is solid, its elastic strength will 
likely be unimportant for the equilibrium tide. However, the density contrast between the 
high-Z core and its surroundings will affect its tidal distortion. For simplicity, consider 
the equilibrium tide in a nonrotating "planet" composed of two incompressible fluids having 
different densities: p = p 2 in the core, < R < R c , and p = p\ < p<i outside it, R c < R < R p . 
As usual, the perturbing tidal potential is quadrupolar, $i >ex t = eR 2 P2(cos9). Since the 
vorticity vanishes except at the interface between the two fluids, the displacements can be 
assumed to be proportional to the gradient of a scalar, £ = V%, where x satisfies Laplace's 
equation and is oc P 2 (cos9); x is discontinuous at the interface but its radial derivative must 
be continuous. This idealized problem can then be worked out analytically by matching the 
radial parts of x an d of $i jSe if across the interface. Included among these conditions is that 
the perturbed interface should remain an equipotential, 

dR~ GMJRI atK - K ^ 

which is analogous to the free boundary condition at the surface. If 5 = pij p\ > 1 is the 
ratio of densities and 77 = R c /R p < 1, it can then be shown that the radial strain in the core 
is related to the radial strain at the interface by 



R 



5 + (5 - lW U 
R<Rc 55-3(5-1X1-^) R 



(50) 
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In the limit 77 -C 1 (recall that we expect R c < 0.2i? p ), the factor on the righthand side 
reduces approximately to 5/(25 + 3). Since 5 ~ 2-4, the core distorts substantially less 
than the main body of the planet, though it does not remain perfectly spherical. We would 
therefore still expect the tide to generate short- wavelength inertial waves in a rotating planet 
with a core, but compared to the estimate (138|) made for a rigid core, the tidal Q would 
increase by a factor ps [1 — 5/(25 + 3)]~ 2 = [(26 + S)/2(S — l)] 2 , i.e. by one half to one order 
of magnitude. 



5. Summary 

Motivated by poss ible app licat i ons to short-period extrasolar planets, and by past work 



by lOgilvie fe Linl (120041 ) and by lWul (l2005bl ). we have studied the dynamical tide in isentropic 
fluid bodies with and without cores. Our goal has not been to obtain precise numerical results 
for realistic planetary structures, since these are still quite uncertain, but rather to provide a 
simple yet semi-quantitative physical explanation of why a rigid core should give rise to short- 
wavelength inertial waves. We do this essentially by a combination of WKB and perturbation 
theory, in which the small parameters are (i) wavelength over radius, and (ii) core radius 
over planetary radius. The essential element in our model is the non-specular reflection of 
inertial waves at a surface. Such reflections can cause dramatic changes in wavelength when 
the surface is nearly perpendicular to one of the directions of the wavevector allowed by 
the WKB dispersion relation at the tidal frequency. In other words, we emphasize critical 
latitudes rather than wave attractors. The spirit of our analysis is very much more local 
and informal than that of most previous work on modes and tides in rotating bodies; we 
hope that the local approach will be accepted as complementary rather than contradictory 
to global analyses. 

In order to obtain an upper bound on the tidal Q from our local approach, we consider 
the production of waves so close to the critical latitudes, and hence so short in wavelength, 
that they damp nonlinearly after a single encounter with the core. Our assumption is that 
wave attractors, which would involve waves that encounter the core and the surface repeat- 
edly before damping, can only lower Q further, as would the escape of inertial waves into 
the stably stratified radiative zone near the surface. We have also examined the physical 
basis for the assumption of a rigid core. We find that rock would probably be molten at the 
central temperatures and pressures expected for hot Jupiters, and that even if the core were 
in a solid phase, it would be sufficiently plastic that its large-scale tidal distortion would 
closely approximate that of a fluid. However, because of the core's self-gravity and higher 
density, it will distort less than its surroundings, so that short-wavelength inertial waves will 
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still be excited. Our upper bound to the tidal Q is in the range 10 7 -10 8 (i? c /0.2i? P 



We have criticized the tidal calculations by IWul (l2005bl ) and by llvanov fc Papaloizou 



(120071 ) for coreless isentropic bodies. We believe that their results, though probably correct 
for bodies that are confined within a rigid outer boundary, overestimate the excitation of 
short inertial waves — and therefore underestimate Q — for the astrophysically relevant case 
of a free outer boundary. We support our case in part by reference to explicit analytical 
solutions of the tidal response in idealized coreless models, albeit ones that are themselves 
not entirely realistic. Since short inertia l waves are certainly p r evalen t a mong the linear 



modes of such bodies, as shown vividly by lRieutord fc Valdettard (119971 ) and lRieutord et al. 



( I200ll ). and since they can be excited in some circumstances (e.g. when a core is present), 
a general theorem concerning the conditions under which short-wavelength modes can and 
cannot be forced by external potentials would be desirable. 
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paper, though we do not suggest that they endorse all of its conclusions. We thank Adam 
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A. Appendix: Tidal Q for an eccentric, synchronous orbit 



Through first order in orbital eccentricity e, the quadrupolar part of the ti dal potenti al 
acting on a rotationally aligned but not necessarily synchronous planet is (e.g. IZahnlll977l ) 



$1 



cxt 



a" 



+e 




-i2(Cl-n)t 



Y 2;0 (9)e- int + ^F 2 ,_ 2 



-i(2Q-n)t _ Jy^Y, 



2,21 



-i(3n-20)t 



• (Al) 



Here the azimuthal coordinate corotates with the planet. Spherical harmonics rather 
than Legendre functions have been used, to clarify the relative strengths of the various 
components. For circular orbits and non-synchronous spins, only the second term in curly 
braces contributes to dissipation, because the rest vanish or are constant in time. This 
appendix is devoted to synchronous (n = fl) but slightly noncircular cases, < e C 1. Then 
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all of the variable parts of the perturbing potential are oc e and have the same frequency: 



$ 1 = e L-MlR*il — Real 




47T. 
5" 



7^ 



2,2^, . 



-int 



(A2) 



where the factor / accounts for the self-gravity of the equilibrium tide as in eq. (|28|) . Fol- 
lowing the prescription of §3.2.2[ the corresponding amplitudes of the incoming component 
of the equilibrum response at the core are 



A- 



(m) 
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(A3) 



m = 0. 



Notice that the m = component is the same, apart from the factor of eccentricity, as that 
of the \m\ =2 component in the nonsynchronous circular case. The radial energy flux of 
outgoing waves near the critical latitudes at the core is 
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The wave power f)33p then becomes 



E(Ad) 
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(A5) 



The sum of the maximum potential energies in the tidal distortions associated with each of 
the variable components in eq. (1A2j) is 



„ KGMlRl 
AE ss 21e 2 — 



(A6) 



Therefore, using the same approximations that led to eq. (|38|) . the tidal Q is 



Q P3 2.1 x 10 7 



P 
Id 



Mi 



Pi 



P r 



0.2/2, 



-5 



0.1 rad 



-1 



(A7) 



This is very close to the result (l3"8i) [but slightly different by virtue of the rational numbers 
entering eqs. (1331) . (|37|) . (1A5[) . and (1A6|) ] because both the synchronization and circularization 
tides are dominated by one of their harmonic components. 
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